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We calculate the neutron electric dipole form factor induced by the CP violating 
#-term of QCD, within a perturbative chiral quark model which includes pion and 
kaon clouds. On this basis we derive the neutron electric dipole moment and the 
I electron- neutron Schiff moment. From the existing experimental upper limits on 

the neutron electric dipole moment we extract constraints on the ^-parameter and 
compare our results with other approaches. 

in 
o 

a: 

CLT I. INTRODUCTION 

The understanding of CP violation is one of the long-standing and challenging problems 
of particle physics. So far the effects of CP violation have only been observed in the K 
and B hadron systems |l| and are in good agreement with the predictions of the Standard 
Model (SM) of electroweak interactions. Among other CP odd observables, a great deal 
of effort has been directed to the study of electric dipole moments (EDM) of leptons, neu- 
trons and neutral atoms. Both SM and various non-SM sources of CP violation have been 



considered (for recent reviews see e.g. Ref. 2], |3|, |4( ). These studies have been particularly 
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stimulated by the expectation of great improvements (2 to 4 orders of magnitude) in the 
experimental sensitivities to EDMs in the next decade (for review see Ref. |3j). 

As it is well known, there are two sources of CP-violation within the SM: the complex 
phase Sckm of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix in the weak 
interaction sector, and the f3-term in the strong interaction sector, which arises due to the 
uon-trivial structure of the QCD vaeuuur Q-fi. The coruple* phase of the CKM uratrix 
provides a consistent explanation of the observed CP odd effects in hadron decays, while 
it gives an imperceptible contribution to the EDMs, far below the sensitivity of present or 
foreseeable experiments. Indeed, the CKM prediction for the neutron EDM, d n , ranges from 
10~ 31 to 10~ 33 e-cm |9] while the present experimental upper limit Q| is 

|4| < 0.63 x KT 25 e-cm. (1) 

On the other hand, CP violation induced by a #-term leads to a sizable electric dipole 
moment for the neutron jlll . ll^ and may significantly contribute to atomic EDMs, while it 
is insignificant for CP violation in hadron decays. The non-observation of the neutron EDM 
and the atomic EDMs imposes a very strict upper bound on the value of 8, of the order 
of 10~ 10 . This unnaturally small value of 6, which is otherwise not restricted by theory, 
is known as the strong CP problem. One elegant solution was proposed by Peccei and 
Quinn Q|, which makes the ^-parameter vanish dynamically. However, the mechanism also 
requires the appearance of a Goldstone boson, the axion, which remains to be discovered. 

Other important contributions to the atomic and neutron EDMs may arise from possible 
physics beyond the SM. In particular, supersymmetric extensions of the SM offer additional 
mechanisms for CP- violation jl^-^^ originating from complex phases in the soft SUSY 
breaking terms and superpotential parameters (for a review see Ref. jl8|). 

In the calculations of atomic and neutron EDMs one faces the problem of translating 
the effect of the CP violation introduced at the quark-gluon level to the processes at the 
hadronic or atomic level. This translation must resort to hadronic and nuclear models as well 
as to a careful treatment of the atomic electron wave functions in the case of atomic EDMs. 
Therefore, confidence of the estimates of the hadronic or atomic CP-violating observables 
depends on reliability of these models. 

The problem of the neutron EDM has been studied within various theoretical approaches: 



current algebra and chiral perturbation theory (ChPT) |l2j, |20j, |2l| , chiral quark models |2 



24j . lattice QCD |25|, |26(, QCD sum rules ]2ly l28l| . an approach based on solutions of 
Schwinger-Dyson and Bethe-Salpeter equations |29j . etc. 

In the present paper we apply to this problem the perturbative chiral quark model 
(PCQM) 1 301 that is a development of chiral quark models with a perturbative treatment 

fin n 

of the pion cloud of nucleon |3J|- 33]. As shown in Ref. [20J, the PCQM is successful in 
description of low-energy properties of light baryons such as the mass spectrum, the electro- 
magnetic, axial and strong form factors, including the quantities which receive a nontrivial 
contribution from the cloud of pseudoscalar mesons: meson-baryon sigma-terms, strangeness 
content of the nucleon, etc. Compared to the models of Refs. |3jJ-|3J| the PCQM contains 
several new features 1 : (i) the SU(3) extension of chiral symmetry in order to include the kaon 
and eta-meson cloud contributions; (ii) consistent formulation of perturbation theory both 
at the quark and baryon levels on the basis of the renormalization techniques and taking 
into account excited quark states in the meson loop diagrams; (iii) incorporated constraints 
from the chiral symmetry (low-energy theorems); (iv) consistency with chiral perturbation 
theory for the case of the chiral expansion of the nucleon mass. 

The purpose of the present work is to calculate within the PCQM not only the neutron 
EDM but also the neutron electric dipole moment form factor (EDFF) induced by the strong 
CP violating ^-term. The neutron EDFF as function of the momentum transfer could be 
the next step in the experimental studies of the CP odd structure of the neutron, after the 
measurement of its EDM. Also, it is known that the atomic EDMs are sensitive to the 
nuclear Schiff moment, which depends on the neutron EDM square radius derived from the 
neutron EDFF. 

This paper is organized as follows. In Sect. [H]we give a brief introduction to the PCQM. 
Sect. IIHI deals with the calculation of the electric dipole form factor, the neutron EDM 
and the electron-neutron Schiff moment induced by strong CP violation. Here, we extract 
the constraints for the CP violating ^-parameter from the existing experimental data on 
the neutron EDM and compare our results to some other approaches. Finally, we give a 
summary of our results and conclusions. 



For details sec Ref. 
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II. THE PERTURB ATIVE CHIRAL QUARK MODEL 



The basis of the perturbative chiral quark model (PCQM) [3Cj is an effective chiral 
Lagrangian describing the valence quarks of baryons as relativistic fermions moving in an 
external field (static potential) V e s(r) = S(r) + j°V(r) with r = \x\ , which in the SU(3)- 
flavor version are supplemented by a cloud of Goldstone bosons (n, K, if). Treating Goldstone 
fields as small fluctuations around the three-quark core, the linearized effective Lagrangian 
is written as: 

C cS (x) = q(x)[i ?-S(r)-j V(r))q(x) + lj2lW*)] 2 

1 i=i 

+ Cf il \x) + C xS B(x). (2) 



Here we defined 



Cf (1 \x) = -g(*)* 7 5 ^W)g(*) . (3) 



The additional term £ x sb in Eq. (0) contains the mass contributions both for quarks and 
mesons, which explicitly break chiral symmetry: 

£ x sb(x) = -q(x)Mq(x) - -Tr[&(x)M] . (4) 

8 

Here, $ = is the octet matrix of pseudoscalar mesons, F = 88 MeV is the pion decay 

i=l 

constant in the chiral limit, A4 = diagjm, rh, m s } is the mass matrix of current quarks (we 
restrict to the isospin symmetry limit m u = rrid = rh) and B = — < 0\uu\0 > / F 2 is the 
quark condensate constant. We rely on the standard picture of chiral symmetry breaking 
and for the masses of pseudoscalar mesons we use the leading term in their chiral expansion 
(i.e. linear in the current quark mass): 

Ml = 2mB, M 2 K = (rh + m s )B, M 2 = -(rh + 2m s )B . (5) 

o 

In our analysis we use the following set of parameters: 

m = 7 MeV, m s = 25m, B = M^ + /2m = 1.4 GeV . (6) 

The meson masses satisfy the Gell-Mann-Oakes-Renner and Gell-Mann-Okubo relations. 
In addition, the linearized effective Lagrangian fulfills PCAC. The properties of baryons, 
which are modeled as bound states of valence quarks surrounded by a meson cloud, are then 
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derived using perturbation theory. At zeroth order, the unperturbed Lagrangian simply 
describes a nucleon as three relativistic valence quarks which are confined by an effective 
one-body static potential V e g(r) in the Dirac equation. We denote the unperturbed three- 
quark ground-state as \4>q), with the normalization (0o|0o) — 1- We expand the quark field 
q in the basis of eigenstates generated by this potential as 

Qi x ) = J2 b a u a (x) exp(-i£ a t) (7) 

a 

where the quark wave functions {u a } in orbits a are the solutions of the Dirac equa- 
tion including the potential V e g(r). The expansion coefficients b a are the corresponding 
single quark annihilation operators. All calculations are performed at an order of accu- 
racy o(l/F 2 , rh, m s ). In the calculation of matrix elements, we project the quark diagrams 
on the respective baryon states. The baryon states are conventionally set up by the product 
of SU(6) spin-flavor and SU(3) C color wave functions, where the nonrelativistic single quark 
spin wave function is replaced by the relativistic solution u a (x) of the Dirac equation. 

In our description of baryons we use an effective potential Ks ( r ) — S(r) + 7 c V(r) which 
is given by a sum of a scalar potential S(r) providing confinement and the time component 
of a vector potential 7°V A (r). Obviously, other possible Lorenz structures (e.g., pseudoscalar 
or axial) are excluded by symmetry principles. It is known from lattice simulations that 
a scalar potential should be a linearly rising one and the vector potential is thought to be 
responsible for short-range fluctuations of the gluon field configurations js^J. In our study 
we approximate V e g(r) by a relativistic harmonic oscillator potential with a quadratic radial 
dependence j3] 

S(r) = M 1 + cir 2 , V(r) = M 2 + c 2 r 2 . (8) 

The model potential defines unperturbed wave functions for the quarks, which are subse- 
quently used to calculate baryon properties. This potential has no direct connection to 
the underlying physical picture and is thought to serve as an approximation of a realistic 
potential. Note, that this type of the potential was extensively used in chiral potential mod- 
els 1311-1331. A positive feature of this potential is that most of the calculations can be done 

fin n 

analytically. As was shown in Refs. |3J|-[3i2] and later on also checked in the PCQM 30], 
this effective potential gives a reasonable description of baryon properties and can be treated 
as a phenomenological approximation of the long-range potential dictated by QCD. 
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The use of a variational Gaussian ansatz for the effective potential (JHJ) gives the following 
solution for the ground state (for the excited quark states we proceed by analogy): 

/ 

Uq\x) = N exp 



2R 2 




Xs Xf Xc , (9) 



where N = [7r 3 / 2 i? 3 (l + 3p 2 /2)] -1//2 is a normalization constant; Xs, Xf> Xc are the spin, flavor 
and color quark wave functions, respectively. The parameter p, setting the strength of the 
"small component", can be related to the axial charge qa of the nucleon. In the leading 
order (3-quark-core) approximation, this relation is [3j| 

The parameters of the effective potential V c s(r) can also be expressed in terms of p and R: 

Here, Sq is the single-quark ground-state energy In our calculations we use the value 
(7^4=1.25. Therefore, we have only one free parameter in the model. In our numerical 
study, R is varied in the region from 0.55 fm to 0.65 fm, which is set and constrained by 



30 



nucleon phenomenology 3Jj | . Such a variation of the parameter R slightly changes the phys- 



. In this paper we also test a sensitivity of the neutron EDM to 



ical quantities up to 5% 
a variation of R. 

The expectation value of an operator A is defined as 

oo -n r „ 

(A) = B (0 O | J2 ~ / d4 xi ■■■ d 4 x n T[Cj( Xl ) . ..dixJAWfa)* (12) 

n=l n - 

where Lj is the full interaction Lagrangian which may contain both CP-even and CP-odd 
terms, as discussed below. The superscript " B" in Eq. (fT2"|) indicates that the matrix ele- 
ments are projected on the respective baryon states and the subscript "c" refers to contri- 
butions from connected graphs only. 

For the evaluation of Eq. (|12|) we apply Wick's theorem with the appropriate propagators 
for quarks and mesons. For the quark field we use a vacuum Feynman propagator for a 
fermion in a binding potential. In the calculation of meson-quark loops we include only the 
ground state in the quark propagator, which leads to the following truncated form: 

iG q (x,y) = (0\T{q(x)q(y)}\0) - 6{x - y^u^u^e-^^ . (13) 
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Note, that in our previous papers we estimated explicitly the contribution of the low-lying 
excited quark states in the quark propagator to the physical quantities. Their contribution is 
about ~ 10 — 15% with respect to the ground state contribution. Therefore, a restriction of 
the quark propagator to the ground states is a reasonable approximation. For completeness 
we also calculate the corrections to the neutron EDM due to the inclusion of excited quark 
states: the first p-states (lpi/2 and lj>3/2 in the non-relativistic notation) and the second 
excited states (1^3/2,1^5/2 an d 2sj/ 2 ), i.e. we restrict to the low-lying excited states with 
energies smaller than the typical scale of A = 1 GeV of low-energy approaches. 
For the meson fields we use their free Feynman propagators: 

iAv{x -y) = (0\T{M^(y)}\0) = S tJ J { ^H^-^-t ] - {U) 

III. STRONG CP VIOLATION AND THE NEUTRON EDM 

In this section we study the electric dipole moment form factor (EDFF) and the electric 
dipole moment (EDM) of the neutron induced by the #-term, using our perturbative chiral 



guar 



2l|- 



£ model (PCQM). We follow a chiral approach for the treatment of the #-term 



24j and then apply the PCQM for the calculation of hadronic matrix elements. 



At the fundamental level, the QCD Lagrangian is: 



Cqcd = -\ tr(GV GT) + qiif- M)q + ^tr{G ^G^) . (15) 



Here, the last term is the CP-violating #-term that cannot be eliminated due to the nontrivial 
topology of the QCD vacuum. As usual, T>^ is the covariant derivative, G^ v is the gluon 
stress tensor (in SU(3) matrix notation) and G^ = {l/2)e^ up G (Tp is its dual tensor. As 
it is well known, doing a chiral U(l) transformation in flavor space, one can remove the 
gluonic ^-term from the Lagrangian (|TK|) and pass it as a (CP-violating) complex phase to 
the quark mass operators. For a small value of 9 the CP violating term becomes (for details 



see Refs. 



in 
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ad); 

£cpv = i0fhq(x)j 5 q(x) (16) 



where q(x) denotes a flavor triplet, and, consequently, this term is a flavor-5'L r (3) singlet. 
The mass coefficient is: 

m u m d m s 

m = , (17) 

m u m d + m u m s + m d m s 
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which would vanish if any of the flavors were massless. From this term we construct the 
effective chiral Lagrangian. As usual, we do this by introducing the chiral field e % ^l F and 
expand it in powers of 

/ $\ _ _ _$ 

£cpv = i0mqi5 exp I ry 5 — J q = iOmq^q - Qmq—q + 0(<3> ) . (18) 

It turns out that the first term in the r.h.s. of Eq. (|T8|) does not contribute to the EDFF at 
one loop because the diagrams involving this vertex do not contain the spin-flip structure 
ctjv • q, the product of the neutron spin operator (Tn and of the 3-momentum of the photon 
q. The leading contribution to the EDFF comes from the linear term of the expansion 
in Eq. ([18)1 . Neglecting higher order terms, we adopt for our analysis the CP- violating 
interaction in the form: 

C s ^ = -6mqjq- (19) 
To guarantee electromagnetic gauge invariance in non-covariant approaches (see discussion 
in Refs. Q Q) 

we have to work in the Breit frame, where the momenta of the initial 
neutron, final neutron and photon are defined respectively as: 

p = (E, -q/2) , p' = (E, q/2) , q = p' - p = (0, q ) . (20) 



Here, E = \lrn? N + (p/4 is the nucleon energy, uin is the nucleon mass, and q 2 = —Q 2 = —q 2 
is the momentum transfer squared. The neutron EDFF, D n (Q 2 ), is defined in the standard 
way through the neutron matrix element of the electromagnetic current: 

(n(p')\ J„(0) |n(p)> = u n (p')U,Fl(Q 2 ) + ^a^q»F 2 (Q 2 ) (21) 



°»vlzq v D n (Q 2 ) + (7 M g 2 - 2m N q^)^A n {Q 2 ) 



Unip) 



where, in addition, F^(Q 2 ) and F 2 (Q 2 ) are the well-known CP-even neutron electromagnetic 
form factors and A n (Q 2 ) is the neutron anapole moment form factor. In our model at one 
loop level the neutron EDFF, D n (Q 2 ), is given in the Breit frame by 

E xl sf i*N ■ qX Ns D n (Q 2 ) = y \^ J 5{t)d A xd A Xl d 4 x 2 e-^ (22) 



m N 

x T 



where xn s and Xn , are ^he nucleon spin wave function (w.f.) in the initial and final state, o?q 
is the nucleon spin operator, jrf^^ anc l Copy are the linearized parts of the CP-even and 
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CP-odd Lagrangian parts [c.f. Eqs. (J3J) and (fTU)). respectively] describing the strong interac- 
tions of pseudoscalar mesons and quarks; Jo is the time-component of the electromagnetic 
current of the charged pseudoscalar meson fields: 



= e (tT id^ + + R- id fi K + ) + h.c. 



(23) 



The only non- vanishing diagram contributing to the strong CP- violating part of the neutron 
EDFF is the so-called meson-cloud diagram shown in Fig.l. This diagram has two contribu- 



tions: from pion and from kaon loops. T 



different theoretical approaches 



re pion 



21. 



23 



cloud effects have already been estimated in 



24j . In addition to the pion cloud the kaon 



loop effects have been calculated in the framework of Heavy Baryon Chiral Perturbation 
Theory HBChPT Here we include both pion and kaon cloud contributions. 

The result for the strong CP-violating contribution to the neutron EDFF is given by 



D 



str:ir 



(Q 2 



D 



str:K 



D 



str:& 



(Q 2 



str 9-kNN 9ixNN 



2E 



(Q 2 ) 

d 3 k F nNN [{k + qf] F vNN (P 



(27T) 1 



w$(k + q) w$(k) w$(k + q) + w<$>{k 



(24) 



#5) 



where $ = 7r or K, w$(q) = y Mf + q 2 is the meson energy, and SU(6) spin-flavor 

factor, which is c^ r = 1 for the pion- loop diagram and cf£ r = 1/5 for the kaon- loop diagram. 
Notice that the expression for the neutron EDFF is written in terms of the strong CP- 
conserving g n NN and CP- violating g^NN pion-nucleon coupling constants, which are basic 
quantities in the effective low-energy pion-nucleon Lagrangian derived from QCD [8j. The 
strong coupling constant g n NN satisfies the Goldberger-Treiman relation: 



m N 

g?rNN — gA -p-, 

where g^ is the axial nucleon charge, whereas the CP-violating coup 
to the result derived in the context of the chiral quark model of Ref. 



(26) 



ing g n NN is identical 
2J: 



gnNN 



m 



(27) 



Here 7 is the isovector-scalar two-quark condensate in the nucleon: 



(iV|gr 3 g|iV) = ju n t 3 u n 



(28) 
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This factor 7 coincides with the so-called relativistic reduction factor j^. In the PCQM, 
7 = 5/8 [30(. Finally, F nNN and F nNN are the normalized strong CP-conserving and CP- 
violating irNN form factors, which regularize the divergent loop integral. In the PCQM 
these form factors are given by 2i3| 

k 2 R 2 



F nNN (k 2 ) = exp(-k 2 R 2 /A) 
and 

K NN (k 2 ) = exp(-k 2 R 2 /4) 



1 + 



1 + 



3&4 

k 2 R 2 1 - 3g A /5 



(29) 



(30) 



4 W5-1. 

where R = 0.6 fm is the PCQM dimensional parameter defining the quark wave function. 
The numerical results for the strong CP-violating contributions from the pion and kaon 
loops to the neutron EDFF as functions of the momentum transfer are displayed in Fig.2. 
Here, the following comment is in order. The lack of covariance in our model makes the 
reliability of the above results for the neutron EDMFF decreasing in the region of large 
momentum transfer, where the relativistic effects become non-negligible. Their typical size 
is determined by the ratio q 2 / (4m 2 N ), where q is the three-momentum transfer. Therefore, 
the relativistic corrections are expected to be linearly growing with Q 2 = q 2 . However, in 
the region Q 2 = q 2 < 0.4 GeV 2 they do not exceed ~ 10% and can be neglected. 
The neutron EDM, d n , is defined as the value of the neutron EDFF at zero recoil: 

7str;$ _ n str ^(r\\ — „ „str 9ttNN9ttNN f d 3 k F n NN(k 2 ) F n NN(k' 



2m N J [2ixy w%(k) 
Now, it is worth checking if our result for the pion-cloud contribution to d^ r ' ,7T in Eq. (j3*Tj) is 



12j for the leading term 



consistent with the model-independent prediction derived in Ref. 
in the chiral expansion: 

- tr;7r = 9,NN ^NN ™V 

AiT 2 m N B M n K ' 

To this end, we drop the normalized form factors F nNN and F nNN in Eq. (|3~Tj) by substitut- 
ing Ftttvat = F n NN = 1, and analyze this equation using alternatively cutoff or dimensional 
regularizations. Both methods of regularization give the same result, which also coincide 
with Eq. (|32p. Therefore, our approach is consistent with QCD in the local limit, when 
F-n-NN — F nNN = 1. The nontrivial form factors F nNN and F wNN provide an ultravio- 
let convergence for the EDM. The leading contribution of the kaon-cloud diagram is also 
proportional to the chiral logarithm but contains a model-dependent coefficient cf£ r = 1/5: 

Jstr;K _ str 9nNN 9nNN , m N /oox 
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We remark that the coefficient c s S r was calculated previously in Heavy Baryon Chiral Per- 

n 

turbation Theory (HBChPT) 21], where it was expressed through the parameters of the 
chiral Lagrangian: 



c 



str = D-Fb F -b D 

K D + F b F + b L 



(34) 

>D 

Here D and F are the axial- vector couplings; bp, and bp are the low-energy constants. Using 
the actual values [2]] of D = 0.80 and F = 0.46 fixed from a fit of semileptonic hyperon 
decays, and 6p = 0.079 GeV -1 , b F = —0.316 GeV -1 determined from the calculation of 
baryon masses and the ttN sigma-term up to fourth order in the chiral expansion, we deduce 
the prediction of HBChPT of cf£ r = 0.45. This is more than a factor two larger than the 
prediction of our model. 

Our results for the neutron EDM induced by the (9-term at one-loop approximation are: 

(a) partial pion and kaon loop contributions 

<ff~ ;7r = 1.37 x 10 -16 x 6 [e ■ cm] , df' K = 0.05 x 10~ 16 x 6 [e • cm] . (35) 

(b) the total pion and kaon loop contribution 

d str = d str;n + = 1M x 1Q -16 x Q [ g . cm ] _ ( 36 ) 

As seen from Eqs. f!35j) . the kaon contribution to the neutron EDM is smaller than the pion 
contribution by a factor ~ 28 . 

Now, let us estimate the sensitivity of our results for the neutron EDM on the model 
approximations: (i) variation of the size parameter R from 0.55 to 0.65 fm, (ii) contribution 
of the excited quark states and (iii) two-loop corrections. 

(i) For the range of values R = 0.6 ± 0.05 fm the values of the quantities, defined in 
Eqs. f!35jl and (JHHJ), are varying within the following ranges: 

df = (1.42 ± 0.1) xl0" 16 x# [e-cm], df ;7r = (1.37 ± 0.1) x 10~ 16 x 9 [e ■ cm](37) 
df' K = (0.05 ± 0.01) x 10~ 16 x 6 [e • cm] . 

(ii) The contribution of the excited quark states to the quark propagator, neglected in 
the present study, according to our estimation is of the order of 10% both for the pion and 
kaon meson cloud diagrams. This result is in agreement with the previous studies of the 
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electromagnetic nucleon form factors, the electromagnetic N — > A7 transition and the axial 
3fl- 



nucleon form factor 

(iii) The two-loop corrections to the one-loop result can be roughly estimated using a 
naive dimensional analysis |2a], from which it follows that their contribution is suppressed 



by a factor 

/ A \ - / 1 



e 



) = {^fr)^ oa (38) 



V4vrF/ KAttFRs 

and, therefore, can be safely neglected at the level of accuracy adopted in the present 
analysis. 

The current experimental bound on the neutron EDM in Eq. JIJ used in Eq. (J3fij) leads 
to the following upper limit for the QCD angle 9: 

\9\ < 4.4 x 1(T 10 . (39) 

The results of our model are in a reasonable agreement with the predictions of other 
chiral approaches shown in Table I. 

It is worth noticing that our prediction for the kaon cloud contribution d^ 7 "^ = 0.05 x 
10~ 16 x 9 [e ■ cm] is much smaller than the analogous result of HBChPT V2W: di tr ' K = 
1.1 x HT 16 x 9 [e-cm]. 

Another CP-odd parameter which may have important implications for CP violation in 
atoms is the electron-neutron Schiff moment S', defined in terms of the neutron EDFF 
D n {Q 2 ) as Q 

~dD n {Q 2 ) 



S' 



dQ 



2 



(40) 



Q 2 =0 



Using our present calculation for D^ r (Q 2 ) induced by strong CP- violation, we have 

nlstr;® s tr 9ttNN9ttNN f d 3 k F^NN (k 2 ) F^NN (k 2 ^ 

O — 6 C/f 



* 2m N J (2vr) 3 ^(fc) 

3 5 k 2 1 M% F> NN (k 2 ) 2 p F^ NN (P 

Awl(k) QwUk) 8m 2 N w%{k) F„ NN {P) 3 F nNN (k- 



(41) 



where F^ NN and F" NN are the first and second derivatives of the F^tvv form factor with 
respect to k 2 . 

Again, as in the case of the neutron EDM, we may check the consistency of our approach 
with ChPT at leading order in the chiral expansion. To this end we drop the normalized 
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form factors F nNN and F nNN in Eq. (J4*T]) substituting F nNN = F nNN = 1 and, therefore, 
F^nn = ^"nn = 0- Then, keeping the leading terms in the chiral expansion (which are 
ultraviolet-convergent) we get: 



qlstr;® _ str 9nNN9nNN 



d 3 k 



5 k' 



2m N J (2tt)3 w l(k) U 6 w |(£). 
Carrying out the integration in Eq. (J42j) explicitly, we arrive at the expression: 



(42) 



olstr;$ _ str 9ttNN 9-kNN 

* A8n 2 m N Ml ' 1 ' 

In the case of the pion-cloud contribution (c** r = 1) this result coincides with the leading 
order result of ChPT 

Here, in the numerical estimate of S" we use the values g 2 NN /4Tr rs 14 and (J^nn = 0.027 x 
from Ref. Q 

Finally, from Eq. (J4*T]) we obtain the following results for the electron-neutron Schiff 
moment S' in our model: 

(a) partial pion and kaon loop contributions 

s ,str;n = 3 g g x 1Q -2 x [ e / GeV 3 ] , S"^;* = 0.05 x 10~ 2 x [e/GeV 3 ] . (45) 

(b) the total pion and kaon loop contribution 

S ,str = S ,str;ir + g> str;K = 3 gl x 1Q -2 x £ [g/QeV 3 ] . (46) 

The above results have been derived in one-loop approximation, with the quark propagator 
truncated to include only the ground state and for the size parameter value R = 0.6 fm. 

As seen from Eq. (I45|) . the kaon contribution to the electron-neutron Schiff moment is 
suppressed by a factor ~ 10~ 2 compared to the pion contribution. Comparing our prediction 
for S' in Eq. (J46|) with the leading order prediction of ChPT j^J, shown in Eq. (J44j) we 
conclude that both results are numerically rather close. 



IV. SUMMARY AND CONCLUSIONS 



In this work we applied the perturbative chiral quark model to the calculation of the neu- 
tron electric dipole form factor induced by a strong CP violating #-term. We have taken into 
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account both pion and kaon cloud contributions. From the existing experimental constraints 
on the neutron electric dipole moment we have derived an upper limit on the CP violating 
parameter 9, which is compatible with the corresponding limits of other existing approaches. 
However, we have found that our prediction for the kaon cloud contribution to the neutron 
EDM is much smaller than the analogous result from Heavy Baryon Chiral Perturbation 
Theory (HBChPT). We have also calculated the electron-neutron Schiff moment, a quantity 
that can be used for the calculations of the electron-nucleus Schiff moments on the basis of 
specific nuclear models. We have found that our result is numerically consistent with the 
leading order prediction from chiral perturbation theory. 

In the present paper we have not considered the sources of the CP-violation beyond the 
SM, which, however, may have an important impact on the neutron electric dipole moment 
form factors. The corresponding results in our model will be published elsewhere. 
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r i 

FIG. 1: The strong CP-violating contribution to the neutron electric dipole moment form factor: 
meson cloud diagram. The vertex denoted by a black filled circle corresponds to the insertion of 
the strong CP violating interaction of Eq. (1191) . 




FIG. 2: The contributions from pion and kaon loops to the neutron EDM form factor induced by 
strong CP violating 6 term. 
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TABLE I: Theoretical estimates of the neutron EDM induced by strong CP violation in units of 
101 x 1(T 16 e-cm. 



K tr \ 


Model 


Reference 


2.7 


MIT bag model 


Baluni [11] 


3.6 


Current algebra 


Crewther et al. [12] 


3.3 


Effective chiral approach 


Pich et al. [20] 


6.7 


HBChPT 


Borasoy [211 


3.0 


Chiral bag model 


Musakhanov et al. [22] 


1.4 


Cloudy bag model 


Morgan et al. [23] 


1.17 


Chiral quark-meson model 


McGovern et al. [24] 


2.4 


QCD sum rules 


Pospelov et al. [28] 


1.42 


This work 





